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RESERVOIR  SYSTEM  ANALYSIS  FOR  WATER  QUALITY 


By 

JAMES  H.  DUKE*/,  DONALD  J.  SMITH^/  AND  R.  G.  WILLEyI^ 
INTRODUCTION 


Water  managers  have  been  attempting  to  operate  reservoirs  In  a  'best 
water  management  strategy"  to  meet  specific  project  objectives  since  the  first 
man-made  impoundment  was  constructed.  Two  things,  however,  require  water 
managers  to  occasionally  reevaluate  a  “best  water  management"  strategy.  As 
time  goes  on,  project  objectives  may  change  due  to  either  an  individual 
owner/operator's  needs  and  desires  or,  in  the  case  of  large  projects,  due  to 
public  Interests  changing  or  finally  being  recognized.  Secondly, 
state-of-the-art  techniques  for  impoundment  regulation  change,  and  may  require 
a  reevaluatlon  of  operational  impacts  due  to  structural  or  nonstructural 
changes . 


STATE-OF-THE-ART 

Wtile  water  managers  are  probably  performing  a  best  project  operation  for 
w. .:•••*  quality  control  for  the  existing  structural  facility,  studies  should  be 
peri ormed  co  evaluate  possible  improvements  that  could  he  provided  bv 
*i ate-of-the-a rt  structural  modifications  and/or  multi-reservoir  system 
operations.  Evaluation  of  the  water  quality  benefits  due  to  possible 
structural  modifications  can  be  performed  with  several  existing 
state-of-the-art  one-dimensional  computer  programs.  ( 14 ,  .’6 , 7 ,  1  6  ; . 

Computer  programs  to  evaluate  the  impact  on  water  quality  due  to  a 
specific  operation  of  a  large  system  ot  reservoirs  are  not  reidilv  aval  able. 
When  it  w.i--  realized,  several  years  ago,  that  the  U.^.  Arsw  ores  of  Engineers 
must  have  the  capability  to  analyze  the  operation  oi  large  mu . t i reserve! r 
systems  for  water  quality,  the  Hydrologic  Engineering  enter,  lii:  ,  was  t  unded 
to  develop  a  computer  program  to  meet  tills  need. 

PROGRAM  DEVELOPMENT 

In  19/8,  various  computer  programs  available  within  the  Cores  tor 
evaluating  reservoir  system  operations  for  water  quantity  were  scr**ece,i  to 
obtain  the  best  generalized  model  tor  adding  water  quality  capability 
[27,12,16).  The  "Simulation  oi  FI  «>,i  tontrnl  and  Conservat  ion  Systems" 
Computer  Program,  HF.C-6,  was  selected  due  to  i-s  generality,  documentation, 
and  level  of  active  support  in  training,  and  maintenance. 


*  Presented  at  the  ASCF.  Hydraulics  liivi  i ■  >n  Specialt”  Conferen  oeur  d'Alene, 

Idaho,  August  14-17,  l‘»84. 

1 J  Consulting  Water  Engineer,  Austin,  Texas,  Member  ASCE 

2J  Resource  Management  Associates,  Lufavette,  California,  Assoc.  Member  ASCE 
3/  Hydraulic  Engineer,  U.S.  Army  Corps  o*  Engineers,  Hvdrologir  Engineering 
Center,  Davis,  California,  Member  ASCE 


The  HEC-5  program  is  designed  to  simulate  the  sequential  operation  ot  a 
reservoir-channel  system  of  a  branching  network  configuration.  Any  time 
interval  from  one  hour  to  a  month  can  be  used.  The  model  contains  the 

capability  to  change  from  one  time  step  to  another  in  order  to  provide  greater 
temporal  resolution  during  certain  periods,  such  as  floods.  Channel  routing 
is  provided  by  any  of  five  hydrologic  routing  techniques.  Reservoirs  operate 
to  (1)  minimize  downstream  flooding;  (2)  evacuate  flood  control  storage  as 
quickly  as  possible;  (3)  provide  for  low  flow  requirements  and  diversions,  and 
(4)  meet  hydropower  requirements.  Hydropower  requirements  can  be  defined  for 
individual  projects  or  for  a  system  of  projects.  Pumped-storage  operation  can 
also  be  simulated.  Sizing  for  conservation  demands  or  storage  can  be 
automatically  performed,  using  the  safe  yield  concept,  and  economic 
computations  can  be  provided  for  hydropower  benefits  and  flood  damage 
evaluation. 

In  1979,  work  was  initiated  to  modify  HEC-5  to  evaluate  reservoir 
operations  for  water  quality  control  in  large  reservoir  systems.  The 
modifications  were  identified  to  be  accomplished  in  three  phases  as  shown  in 
Figure  1.  Phase  I  added  the  capability  to  HEC-5  to  control  water  temperature 
releases  at  one  reservoir  to  provide  for  the  best  combination  of  downstre  in 
needs  at  up  to  three  control  points  (i.e,  river  locations  used  for  controlling 
flow  routing  computations,  and  for  controlling  quantity  and  quality  target 
computations).  The  control  of  water  temperature  is  accomplished  through 
multilevel  intake  structure  operation.  In  September  1979,  a  single  reservoir 
water  temperature  control  program  called  HEC-5Q  [17J  was  completed. 

In  1980,  work  was  initiated  to  modify  the  phase  I  model  to  add  seven  more 
water  quality  parameters  and  capability  to  evaluate  either  two  tandem  (i.e., 
in  series)  or  two  parallel  (i.e.,  on  two  independent  tributaries)  reservoirs. 
In  September  1980,  a  two-reservoir  model  [18]  capable  of  system  operation  for 
three  conservative  and  three  nonconservative  water  quality  parameters,  in 
addition  to  dissolved  oxygen  and  water  temperature  was  completed. 

Following  the  phase  II  development,  the  model  was  modified  with  some 
small  but  significant  additions  and  revisions  during  1981 .  These 
modifications  included  flow  augmentation,  improved  model  efficiency,  and 
recently  developed  selective  withdrawal  routines. 

The  third  and  last  major  phase  of  development  involves  Increasing  the 
HEC-5Q  capability  to  include  up  to  ten  reservoirs  and  thirty  control  points. 
The  HEC  is  currently  testing  the  phase  III  model  on  a  practical  application. 

MODEL  CONCEPTS 


Flow  Simulation  Module 


The  flow  simulation  module  was  developed  to  assist  in  planning  studies 
for  evaluating  proposed  reservoirs  in  a  system  and  to  assist  in  sizing  the 
flood  control  and  conservation  storage  requirements  for  each  project 
recommended  for  the  system.  The  program  can  be  used  In  studies  made 
Immediately  after  the  occurrence  of  a  flood  to  evaluate  preproJecL  conditions 
and  to  show  the  effects  of  existing  and/or  proposed  reservoirs  on  flows  and 
damages  in  the  system.  The  program  should  also  be  useful  in  selecting  the 
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proper  reservoir  releases  throughout  the  system  during  flood  emergencies  m 
order  to  minimi2e  flooding  as  much  as  possible  and  yet  empty  the  system  as 
quickly  as  possible  while  maintaining  a  balance  of  flood  control  storage  among 
the  reservoirs. 

The  above  purposes  are  accomplished  by  simulating  the  sequential 
operation  of  a  system  of  reservoirs  of  any  configuration  for  short  interval 
historical  or  synthetic  floods  or  for  long  duration  nonflood  periods  or  for 
combinations  of  the  two.  Specifically  the  program  may  be  used  to  determine: 

a.  Flood  control  and  conservation  storage  requirements  for  each 
reservoir  in  the  system. 

b.  The  influence  of  a  system  of  reservoirs  on  the  spatial  and  temporal 
distribution  of  runoff  in  a  basin. 

c.  The  evaluation  of  operational  criteria  for  both  flood  control  and 
conservation  (including  hydropower)  for  a  system  of  reservoirs. 

d.  The  expected  annual  flood  damages,  system  costs,  and  system  net 
benefits  for  flood  damage  reduction. 

e.  The  system  of  existing  and  proposed  reservoirs  or  other  alternatives 
including  nonstructural  alternatives  that  result  in  the  maximum  net  flood 
control  benefits  for  the  system  by  making  simulation  runs  for  selected 
alternative  systems. 

Since  many  papers  [19,4,8,9,7,11,23,10]  have  already  been  written 
regarding  the  detailed  technical  concepts  of  the  flow  simulation  module,  the 
remainder  of  this  paper  will  emphasize  the  concepts  of  the  water  quality 
simulation  module. 

Water  Quality  Simulation  Module 

The  water  quality  simulation  module  was  developed  so  that  temperature, 
three  user  selected  conservative  and  three  user  selected  non-conservative 
constituents  can  be  simulated.  The  model  allows  dissolved  oxygen  to  also  be 
simulated  if  the  user  selects  either  carbonaceous  or  nitrogeneous  oxygen 
demanding  constituents,  or  both,  as  one  or  two  of  the  non-conservative 
constituents. 

The  water  quality  simulation  module  accepts  system  flows  generated  by  the 
flow  simulation  module  and  computes  the  distribution  of  all  the  water  quality 
constituents  in  up  to  ten  reservoirs  and  their  associated  downstream  reaches. 
The  ten  reservoirs  may  be  in  any  arbitrary  parallel  and  tandem  configuration. 

The  water  quality  simulation  module  also  selects  the  gate  openings  for 
reservoir  selective  withdrawal  structures  to  meet  user-specified  water  quality 
objectives  at  downstream  control  points.  If  the  objectives  cannot  be 
satisfied,  the  model  will  compute  the  increase  in  flow  (if  any)  necessary  to 
satisfy  the  downstream  objectives.  With  these  capabilities,  the  planner  may 
evaluate  the  effects  on  water  quality  of  proposed  reservoir-stream  system 
modifications  and  determine  how  a  reservoir  intake  structure  should  be 
operated  to  achieve  desired  water  quality  objectives  within  the  system. 
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Each  reservoir  is  assumed  Co  be  a  control  point ,  in  keeping  with  the 
concepts  used  in  the  development  of  the  flow  simulation  module.  The  water 
quality  module  will  allow  for  up  to  thirty  control  points,  including  the 
reservoir  control  points.  The  additional  control  points  may  be  placed  in  the 
stream  system  below  the  reservoirs  at  any  desirable  locations  provided  the 
following  guidelines  are  followed: 

a.  The  most  downstream  point  in  all  systems  must  be  a  control  point. 

b.  The  confluence  of  the  two  streams,  on  which  parallel  reservoirs  are 
located,  must  be  a  control  point. 

c.  The  end  of  the  stream  reach  below  the  more  upstream  reservoir  of  a 
tandem  reservoir  system  (at  the  upstream  end  of  the  more  downstream  reservoir) 
must  be  a  control  point. 

The  water  quality  simulation  module  currently  uses  flow  data  from  the 
flow  simulation  module  at  intervals  of  one  day  and  uses  computational  time 
steps  of  one  day.  Shorter  simulation  time  steps  can  be  used  but  the  model  has 
not  been  tested  with  the  shorter  steps.  The  model  is  limited  to  simulations 
of  one  calendar  year. 

The  reservoirs  are  represented  conceptually  by  series  of  one  dimensional 
horizontal  slices  such  as  those  shown  in  Figure  2.  Each  horizontal  slice  or 
layered  volume  element  is  characterized  by  an  area,  thickness  and  volume.  In 
the  aggregate  the  assemblage  of  layered  volume  elements  is  a  geometric 
representation  in  discretized  form  of  the  prototype  reservoir.  This  one 
dimensional  representation  has  been  shown  to  adequately  represent  water 
quality  conditions  in  many  deep,  well  stratified  reservoirs  by  Eiker  [2],  Baca 
[2]  and  WRE  [28,29,30]. 

Within  each  element,  the  water  is  assumed  to  be  fully  mixed.  This 
implies  that  only  the  vertical  dimension  is  retained  during  the  computation. 
Each  horizontal  layer  is  assumed  to  be  completely  homogeneous  with  all 
isopleths  parallel  to  the  water  surface  both  laterally  and  longitudinally. 
External  inflows  and  withdrawls  occur  as  sources  or  sinks  within  each  layer 
and  are  instantaneously  dispersed  and  homogeneously  mixed  throughout  each 
element  from  the  headwaters  of  the  impoundment  to  the  dam.  It  is  not 
possible,  therefore,  to  look  at  longitudinal  variations  in  water  quality 
constituents.  Module  results  are  most  representative  of  conditions  in  the 
main  reservoir  body. 

Vertical  advection  is  governed  by  the  location  of  inflow  to,  and  outflow 
from,  the  reservoir.  Thus  the  computation  of  the  zones  of  distribution  and 
withdrawal  for  inflows  and  outflows  are  of  considerable  significance  in 
operation  of  the  model.  The  WES  withdrawal  method  [3]  is  used  for  determining 
the  allocation  of  outflow.  The  Debler  inflow  allocation  method  [5]  is  used 
for  the  placement  of  inflows. 

Vertical  advection  is  the  net  interelement  flow  and  is  one  of  two 
transport  mechanisms  used  in  the  module  to  transport  water  quality 
constituents  between  elements.  The  vertical  advection  is  defined  as  the 
interelement  flows  which  result  in  a  continuity  of  flow  in  all  elements. 
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Evaporation 


Outflow 


FIGURE  2 

Geometric  Representation  of  a  Stratified  Reservoir 
and  Mass  Transport  Mechanisms 


Che  other  transport  mechanism  used  in 
constituents  between  elements.  The 


Effective  diffusion  is 
to  transport  water  quality 
diffusion  is  composed  of  molecular  and  turbulent  diffusion  and 
mixing. 


the  module 
effective 
convective 


Wind  and  f  Low  induced  turbulent  diffusion  and  convective  mixing  are  the 
dominant  components  of  effective  diffusion  in  the  epilimnion  of  most 
reservoirs.  In  quiescent  well  stratified  reservoirs,  molecular  diffusion  may 
be  a  significant  component  in  the  metalimnion  and  hypolimnion.  For  deep,  well 
stratified  reservoirs  with  significant  inflows  to  or  withdrawals  from  the 
hypolimnion,  flow  induced  turbulence  in  the  hypolimnion  dominates.  For  weakly 
stratified  reservoirs,  wind  induced  or  wind  and  flow  induced  turbulent 
diffusion  will  be  the  dominant  component  of  the  effective  diffusion  throughout 
the  reservoir.  One  of  two  methods  may  be  selected  by  the  user  to  calculate 
effective  diffusion  coefficients.  For  shallow  weakly  stratified  reservoirs, 
the  wind  controlled  mixing  [15]  method  is  appropriate,  while  the  stability 
method  [15]  is  more  appropriate  for  deeper  well  stratified  reservoirs.  Both 
of  these  methods  have  been  shown  in  numerous  applications  to  adequately 
represent  the  mixing  phenomena  for  heat  and  dissolved  water  quality 
constituents  when  properly  applied. 


The  stream  system  is  represented  conceptually  as  a  linear  network  of 
segments  or  volume  elements.  Each  element  is  characterized  by  length,  width, 
cross-sectional  area,  hydraulic  radius.  Manning’s  n  and  a  flow  and  depth 
relationship  (see  Figure  3).  Flow  rates  at  stream  control  points  are 
calculated  within  the  flow  simulation  module  using  any  one  of  the  several 
programmed  hydrologic  routing  methods.  Within  the  flow  simulation  module, 
incremental  local  flews  (i.e.,  inflow  between  adjacent  control  points)  are 
assumed  deposited  at  the  control  point. 

Within  the  water  quantity  simulation  module,  the  incremental  local  flow 
may  be  divided  into  components  and  placed  at  different  locations  within  the 
stream  reach  (i.e.,  that  portion  of  the  stream  bounded  by  the  two  control 
points).  A  flow  balance  is  used  to  determine  the  flow  rate  at  element 
boundaries.  Any  flow  imbalance  (i.e.,  the  difference  in  the  flow  at  the 
upstream  control  point  plus  all  tributary  inflows  and  the  flow  at  the 
downstream  control  point)  is  distributed  uniformly  to  the  flows  at  each 
element  boundary.  Once  interelement  flows  are  established,  the  depth,  surface 
width,  and  cross  section  area  are  computed  at  each  element  boundary  assuming 
normal  flow. 


Both  the  streams  and  reservoirs  are  represented  by  a  one-dimensional 
assemblage  of  fluid  elements  linked  together  by  interelement  flow  and 
diffusion  (stream  diffusion  is  assumed  to  be  small).  The  interelement  mass 
transports  and  the  fundamental  principle  of  conservation  ot  heat  can  be 
represented  by  the  following  differential  equation  model  of  the  dynamics  of 
temperature  within  each  fluid  element. 
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FIGURE  3 

Geometric  Representation  of  Stream  System  and 
Mass  Transport  Mechanisms 
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T  -  temperature  In  °C 
V  -  volume  of  the  fluid  element  In 
t  -  time  In  seconds 

z  =  space  coordinate  in  meters  (vertical  for  the  reservoir  and 
horizontal  for  the  stream) 

Qz  =*  interelement  flow  in  mVs 

=  element  surface  area  normal  to  the  direction  of  flow  in 
m- 

Dz  =  effective  diffusion  coefficient  in  m?/s 
Q ^  =  lateral  inflow  in  mVs 

Tj  =  inflow  water  temperature  in  degrees  Celsius 
Q0  =  lateral  outflow  in  tn-Vs 
A^  =  element  surface  area  in  m? 

11  =»  external  heat  sources  and  sinks  in  kcal/m~/s 
=  water  density  in  kg/m^ 
c  -  specific  heat  of  water  in  kcal/kg/°C 


This  equation  represents  the  heat  conservation  equation  for  the  fluid 
element.  The  set  of  equations  for  all  elemenfs  within  the  reservoir  or 
stream  system  represents  the  heat  conservation  within  that  system.  All  of 
Die  terms  on  the  right  side  of  Equation  (1)  represent  physical  heat  transfers 
fncluding  the  external  heat  sources  and  sinks.  The  external  heat  sources  and 
sinks  that  are  considered  in  the  module  are  assumed  to  occur  at  the  air-water 
interface.  The  rate  of  heat  transfer  per  unit  of  surface  area  can  be 
expressed  as  the  sum  of  the  following  heat  exchange  components. 


H 

n 


I!  -  H  +  H  -  H  +  H 
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where 


Hn  “net  heat  transfer 

Hs  -=  short  wive  solar  radiation  arriving  at  the  water  surface 
Hsr  -  reflected  short  wave  radiation 
H3  =  long  wave  atmospheric  radiation 
Har  =  reflected  long  wave  radiation 
Hc  =  heat  transfer  due  to  conduction 
Hbr  =  radiation  irom  the  w.iter  surface 
He  =  heat  loss  due  to  evaporation 


All  units  are  In  kcal/m?/s. 


Complete  discussions  of  the  individual  terms  have  been  presented  by 

Anderson  [1|  and  in  Tennessee  Valley  Authority  Report  No.  14  [25J. 

The  method  used  in  the  module  to  evaluate  the  net  rate  of  heat  transfer 
at  the  air-water  interface  was  developed  by  Edinger  and  Geyer  [6].  Their 

method  utilized  the  concepts  of  equiliorium  temperature  and  the  coefficient 
of  surface  heat  exchange.  The  equilibrium  temperature  is  defined  as  the 

water  temperature  at  which  the  net  rate  of  heat  exchange  between  a  water 
surface  and  the  atmosphere  is  zero.  The  coefficient  of  surface  heat  exchange 
is  the  rate  at  which  the  heat  transfer  process  proceeds.  The  equation 

describing  this  relationship  is: 
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(3) 


H  =  K  (T  -  T  ) 
nee  s 

where  Hn  =  net  rate  of  heat  transfer  In  kcal/m2/s 

Ke  =  coefficient  of  surface  heat  exchange  in  kcal/m2/s/°C 
Te  =*  equilibrium  temperature  in  °C 
Ts  =  surface  temperature  in  °C 

A  Heat  Exchange  Program  which  computes  these  terms  is  available  at  the 
HEC  [26]. 

All  heat  transfer  mechanisms  except  short  wave  solar  radiation  apply  at 
the  water  surface.  Short  wave  radiation  penetrates  the  water  surface  and  may 
affect  water  temperatures  several  meters  below  the  surface.  The  depth  of 
penetration  is  a  function  of  adsorption  and  scattering  properties  of  the 
water  [13].  This  phenomenon  is  unimportant  in  the  stream  routines  since 
elements  are  assumed  vertically  mixed. 

In  the  reservoir  routines,  however,  the  short  wave  solar  radiation  may 
penetrate  several  elements.  The  amount  of  heat  which  reaches  each  element  is 
determined  by: 


I  =  (1  -  ?>)  I  e~kz 
o 


(4) 


where 


I  =  light  energy  at  any  depth  in  kcal/m2/s 

S  »  fraction  of  the  radiation  absorbed  in  the  top  font  of  depth 
I0  =  light  energy  at  the  water  surface  in  kcal/iWs 
k  =*  light  extinction  coefficient  in  1/m 
z  =  depth  in  meters 


Combining  equations  (3)  and  (4)  for  the  reservoir  surface  element,  the 
external  heat  source  and  sink  term  becomes: 


H  =  Ke(Te  -  Ts)  -  (1  -  6)  IQe 


-kAz 


(3) 


and  the  external  heat  source  for  all  remaining  reservoir  elements  becomes: 

I  =  I  (1  -  e"kAz) 
z 


(6) 


where 


Iz  =  light  intensity  at  the  top  of  the  element  in  kcal/m2/s 


Water  quality  constituents  other  than  temperature  are  represented  by 
Equation  (1)  with  minor  modifications: 


a.  The  definition  of  the  variable  T  is  generalized  to  represent  the 
concentration  of  any  water  quality  constituent. 

b.  The  distributed  heat  gain/loss  term  A{,H/(oc)  is  : 

(1)  Eliminated  for  conservative  constituents 
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(2)  Replaced  by  a  first  order  kinetic  decay  formulation,  -K^T, 

for  non-conservative  constituents  where  Kj  is  the  decay  rate  in  1/day. 

(3)  Replaced  by  a  first  order  reaeratlon  term,  K.2(D0*-D0),  for 

dissolved  oxygen  where  K.2  is  the  reaeratlon  rate,  DO*  is  the  dissolved 
oxygen  saturation  concentration  at  the  ambient  temperature  and  DO  is  the 
existing  dissolved  oxygen  concentration. 

The  reservoir  reaeration  rate  is  computed  as  follows: 

K  (a  +  bW?)/z  (7) 

where  ~  reaeration  rate  in  1/day  at  20°C 

a,  b  *  empirical  coefficients  derived  by  curve 
fit  from  Kanwisher  [20]  to  be  0.641 
and  0.128  respectively. 

W  »  wind  speed  in  m/s 

z  »  surface  element  thickness  in  meters 

The  stream  reaeratlon  rate  is  computed  using  the  0 ' Connor-Dobbins  [24] 
method : 

K„  =  (D  U)°-V-5  (8) 

?  m 


where  K2  »  reaeratlon  rate  in  1/day  at  20°  C, 

Dm  =  molecular  diffusion  coefficient  in  m^/day 
U  =  flow  velocity  in  m/s 
D  »  average  stream  depth  in  meters 


Ail  first  order  kinetic  rates  are  adjusted 
temperatures  using  a  multiplicative  correction  factor. 

p  =  T  (T-20) 


local  ambient 


where  0  =  kinetic  rate  multiplicative  correction  factor 

Tc  »  empirically  determined  temperature  correction  factor 
T  =  local  ambient  water  temperature 


Reservoir  Solution  Technique 


Within  the  reservoirs,  a  Gaussian  reduction  scheme  is  used  for  solving 
the  differential  equations  which  represent  the  response  of  the  water  quality 
constituents.  Equation  (1)  is  rewritten  in  a  form  where  a  "forward  time  and 
central  difference"  scheme  is  used  to  describe  all  the  derivative  processes. 
For  element  i  adjacent  to  elements  l-l  and  i+1  (see  Figure  4)  the  general 
mass  balance  equation  becomes: 


ELEMENT  i  +  l 


ELEMENT  i-l 


Elements 


V 


■  ho  {[A-/], +0"l}  -  Ti { [-ftF ■] t 


A  I)  ] 

"frj..  +Qd’- 


ji+i 


.!>. ;  ]  .. 
on  •  iv  +  V  +  i.  ..I  I  ~~  I  +od...  r  +  "o  ;  + 

i  +  !  't  I  l  +  l  (.  L  Az  J  1  +  1  x  x 


(10) 


where  subscripts  1,  i-1,  1+1  denote  element  numbers,  and 


V  *  volume  of  the  fluid  element  In 
T  *  I i>-,,iera t ure  In  °C  or  water  quality 
constituent  concentration,  mg/1 
t  -  ■  imputation  time  step  in  seconds 

Az  *  clement  a-ea  at  the  fluid  element  boundary  in  ta- 
D?  *  effective  diffusion  coefficient  in  m^/s 
■z  *  element  thickness  (length  of  stream)  in  meters 
du  ”  upward  advective  flow  (stream  flow)  between  elements  In 

Qj  *  downward  advective  flow  between  elements  in  m^/s 
Qw  *  rate  of  flow  removal  from  the  element  In  m-Vs 
c)x  »  rati'  of  Inflow  to  the  element  in  mVs 
Tx  =  lntlow  watet  temperature  in  °C  or  constituent 
concentration  In  m-./l 

H  ■  external  sources  u.i  sinks  of  heat  in  kcal/s 
=  water  denslcv  in  kc/m^ 
c  =  specific  h--  it  or  water  In  kcal/kg/°C 


Recall  that  the  :/(oc)  terra  is  replaced  by  -K^T  or  Ko  ( uO*-..i' 1 )  for  non¬ 
conservative  water  quality  constituents  and  dissolved  oxygen  respect ivelv .  A 
finite  difference  equation  of  this  type  is  formed  for  each  element  and 
integrated  with  respect  to  time.  The  system  of  finite  difference  mass  balance 
equations  represents  the  response  ot  wite;  quality  within  the  entire  reservoir 
and,  with  t ae  aid  of  numercial  integration  technique,  the  equations  are  solved 
with  respect  to  time. 


Stream  Solution  Technique 


For  the  stream,  t  linear  pro/i  mining  aigo^thm  is  used  iu  solve  i  iullv 
implicit  backward  difference  in  space,  forward  difference  in  rime,  finite 
difference  approximation  of  Equation  (1).  This  approximation  has  the  general 
form 


.t  +  L  .i  +  t+1 

a.  .  ,C.  +  a.  +  a.  .  .C.,  =  b. 

i,i—l  i—l  1,11  i,i+1  i+l  i 


(11) 


where  the  "a"  terms  are  loetli-lents  formed  from  the  area,  dispersion 
coefficient-.,  flows,  lem/hs  of  the  computational  elements  and  time  step  for 
each  volume  element;  the  "C”  terms  ire  the  unknown  temperatures  and  constituent 
concentrations  in  each  volume  element;  the  ”b"  terms  are  constants  formed  from 
initial  conditions  or  previously  computed  conditions,  tributary  inputs  of  heat 
or  mass  loa  is  and,  depending  upon  the  conte <t,  the  reservoir  releases. 
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Two  matrix  formats  are  used  in  the  stream  water  quality  simulation 
module.  The  first  is  used  to  solve  for  temperature  and  constituent 
concentrations  given  all  external  inputs  this  format  is 

! A | c  =  b  (121 

where  |A|  is  the  matrix  of  coefficients,  c  is  the  vector  ;>l  unknown 
temperatures  or  constituent  concentrations  and  b  is  the  vector  of  constants. 
This  first  format  is  used  in  the  water  quality  simulation  module  to  compute 
the  final  results  after  all  reservoir  operations  have  been  completed.  In 
effect,  the  linear  programming  algorithm  is  used  simply  as  a  matrix  solver 
for  a  simulation  model. 


The  second,  and  more  complex,  matrix  format  used  in  the  water  quality 
simulation  module  is  for  determining  the  temperature  and  constituent 
concentrations  that  must  come  from  the  reservoirs  to  satisfy  all  water 
quality  targets  in  the  stream  system.  In  effect,  the  second  format  is  used 
to  (1)  determine  which  control  point  controls  the  release  for  each 
constituent  and  (2)  determine  the  reservoir  release  water  quality  t.  at  most 
closely  satisfies  the  targets  at  the  controlling  point.  This  decision  making 
capability  is  achieved  by  (1)  transforming  the  constituent  concentrations  at 
each  control  point  into  a  specification  of  the  target  and  the  deviation  of 
the  simulated  concentration  above  or  below  the  target,  and  (2)  making  the 
concentrations  in  the  reservoir  releases  unknown  so  that  they  can  be  computed. 


The  transformation  used  at  control  points  to  specify  the  target  is: 


where 


ct+1  =  c  +  ct+1-  ct+1 

l  o  +  L  -i 


ct+1 

1 

s 

n'+l 

+i 

_t+l 


simulated  temperature  in  degrees  Celsius  or  constituent 
concentration  in  mil ligrams/l iter 
target  temperature  in  degrees  Celsius  or  constituent 
concentration  in  milligrams/liter  for  a  control  point 
deviation  of  simulated  temperature  or  constituent 
concentrations  above  the  control  point  target 
deviation  of  simulated  temperature  or  constituent 
concentration  below  the  control  point  target 


(13) 


This  transformation  is  substituted  into  Equation  (11)  to  yield  the  following 
equation  which  is  applied  to  those  volume  elements  that  are  located  at 
control  points. 


„t+l 


(14) 


where  the  (ai,iC0i)  terra  has  been  moved  to  the  right  hand  side  of  the 
equation  since  it  is  known.  Thus,  the  m  x  m  simulation  matrix  has  now  been 
transposed  into  a  m  x  n  rectangular  matrix  where  n  "  ra  ♦  NCP  and  NCP  is  the 
number  of  control  points. 
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Equation  (14)  is  the  general  form  of  the  equation  used  !ur  all  volume 
elements  in  formulating  decision  making  problem.  It  includes,  as  variables, 
the  constituent  concentrations  in  the  reservoir  releases,  although  the 

inclusion  is  not  obvious.  For  those  volume  elments  that  are  just  below 

reservoirs,  the  concentrations  represent  the  constituent  concentrations 

in  the  reservoir  releases.  In  the  simulation  model,  where  the  reservoir 

release  constituent  concentrations  are  known,  the  terms  were 

included  in  the  b  vector  for  those  volume  elements  just  below  reservoirs. 
For  the  decision  making  model,  the  (^C^jterms  are  included  as  unknowns. 
Thus  the  m  x  n  simulation  matrix  has  been  made  even  more  elongated  in 

variables  and  n  is  now  m  +  NCP  +  NRES  where  NRES  is  the  number  of  reservoirs 
in  the  system. 


One  additional  set  of  equations  is  included  in  the  water  quality 
simulation  module  to  ensure  that  realistic  results  are  obtained  in  computing 
reservoir  release  water  quality.  These  equations  are  applied  to  define  the 
range  of  constituent  concentrations  that  may  be  released  from  the 
reservoirs.  Normally  the  range  is  defined  by  two  inequalities: 


> 


rt+l 

min 


(id) 


where 


c 


tH 
;nin 
„t  +  l 
'max 

_t+l 


,  ct+1 

—  max 


(  16) 


minimum  temperature  or  constituent  concentration  in  reservoir 
water  quality  profile 

maximum  temperature  or  constituent  concentra: ion  in  reservoir 
water  quaLity  profile 

final  computed  temperature  or  constituent  concentration  in 
reservoir  release 


In  practice,  these  inequalities  are  written  as  equalities  by  adding  slack  and 
surplus  variables. 


,t+l 
f , 

r 


X 


surplus 


,t+ ! 

i 

min 


(17) 


r:t+1  +  x  ,  =  ct+1 

r  slacK  max 


( 18) 


With  the  problem  so  formulateu,  the jA  matrix  of  Equation  (12)  consists  of 
(m  +  2  *  NKES )  rows  and  (m  +  NCP  +  NKES)  unknowns  and  the  b  vector  consists  of 
(m  +  2  *  NRES)  constants.  ThelAI  matrix  may  be  conceptually  partitioned  as 
shown  in  Figure  5,  where  it  is  assumed  that  reservoirs  are  above  volume 
elements  1  and  3,  that  these  reservoirs  are  in  tandem  and  that  volume  elements 
1,  3,  7  and  m  are  control  points. 
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Stream  Water  Constituent  Concentration  Variables*  Variables  for  Negative  Reservoir  Slack  and 

Deviations  at  Control  Release  Surplus 

_  Points  l  Variables  !  Variables 


i 
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F i GbRL  5 

ructure  of  A,  Matrix  for  Two  Tandem  Reservoir ,  Four  Control  Point  System  With  m  Volume  Elements 


There  are  a  number  of  solutions  that  will  satisfy  a  matrix  that  Is  not 
square  (i.e.,  m  x  .n) .  The  purpose  of  using  a  linear  programming  solution  Is 
to  select  the  solution  that  best  satisfies  the  object!  ves  of  the  reservoir 

operation  on  downstream  water  quality.  However,  it  is  know  that  one  class  of 

solutions  will  never  appear:  at  no  tine  will  the  variables  that  describe  the 
positive  and  the  negative  deviations  from  the  control  point  target 

constituent  concentrat ions  appear  simultaneously  in  the  solution.  At  all 
times,  one  or  the  other  deviation  will  appear  but  not  both.  It  is  also  known 
that  the  reservoir  release  constituent  concentrations  will  always  appear  In 
the  final  solution.  Thus,  selecting  the  solution  that  best  satisfies  the 
objectives  of  the  reservoir  operation  on  downstream  water  quality  reduces  to 
selecting  which  control  point  deviation  variable  appears  in  the  final 
solution  and  the  numerical  value  attached  to  that  variable.  Once  this 
numerical  value  is  known,  It  is  known  that  the  deviation  of  the  opposite  sign 
is  zero  so  that  the  actual  control  point  constituent  concentration  can  be 
computed  using  Equation  (13). 

The  objective  function  Is  used  in  a  linear  programming  formulation  to 
quantitatively  describe  the  desirability  of  any  given  solution  to  a 

formulated  problem.  In  the  water  quality  simulation  module,  a  minimization 
routine  is  used  with  the  objective 

minimize  z  =  p  <:  (19) 

The  actual  value  of  z  Is  Immaterial  to  the  water  quality  simulation  module; 
U  i-  just  an  index  by  which  the  desirability  of  the  solution  is  determined. 
The  vector  c  is  the  same  vector  c  as  in  Equation  (12)  except  that,  as 
previously  described,  it  includes  the  variables  representing 

1.  Hie  deviations  from  the  control  point  targets  for  those  volume 
elements  that  represent  control  points, 

2.  The  constituent  concentrations  In  ail  other  volume  elements,  and 

1.  The  constituent  concentrations  in  the  reservoir  releases. 

The  vector  p  represents  the  penalty  associated  with  the  jppearance_of  a 
given  variable  in  the  final  solution.  Logically,  the  penalties  in  p  are 
nonzero  onLy  at  control  points  and  are  applied  only  for  the  variables  that 
represent  deviations  from  the  target. 

The  water  quality  simulation  module  Is  structured  flexibly  so  that 
different  penalties  can  be  assigned  for  each  control  point,  for  each 
con.stlt.uent  and  for  each  deviation,  above  and  below.  The  magnitude  of  the 
penalty  is  unimportant,  as  long  as  It  Is  nonzero  where  necessary  and 
realistically  represents  the  desired  policy.  For  instance,  for  a  temperature 
target  expressed  as  "the  temperature  at  control  point  1  shall  not  exceed 
20°  C',  or 

T  ;  :>(>  1 20) 

the  penally  for  the  positive  deviation  at  control  point  I  could  be  set  to  1.0 
l ml  (tie  penalty  for  ihe  negative  deviation  could  be  set  to  ;i,u.  Obviously, 
when  trying  to  minimize  z,  Equation  (  1  ^ )  ,  the  linear  programming  algorithm 
would  try  to  ensure  that  the  variable  representing  the  negative  deviation 
would  appear  in  the  final  solution  since  a  lower  valiu  of  the  index  z  would 
result . 


If  it  was  twice  as  important  that  the  temperature  target  at  control 
point  I  be  maintained  chan  at  another  control,  point,  say  J,  then  tne  penalty 
associated  with  a  positive  deviation  from  the  target  at  I  could  be  set  to  2.0 
and  the  penalty  associated  with  a  similar  positive  deviation  at  J  could  oe 
set  to  1.0.  Of  course,  the  penalties  associated  with  negative  deviations  at 
both  I  and  J  would  be  set  to  0.0. 

Similar  logic  is  used  for  setting  penalties  for  constituents  that  must 
always  exceed  a  target  value,  such  as  dissolved  oxygen.  The  nonzero 
penalties  are  applied  to  the  variables  representing  negative  deviations,  and 
the  variables  that  represent  positive  deviations  are  given  penal- its  of  0. 

Gate  Selection 

Once  the  desired  reservoir  release  water  quality  lias  been  computed, 
using  the  above  procedure,  the  water  quality  simulation  module  proceeds  to 
determine:  (1)  the  reservoir  gates  from  which  releases  can  b<_  made;  (2)  the 
gates  that  should  be  used,  and,  (3)  the  water  quality  of  the  releases. 

The  port  selection  algorithm  serves  to  determine  which  ports  should  be 
open  and  what  flow  rate  should  pass  through  each  open  port  in  order  to 
maximize  a  function  of  the  downstream  water  quality  concentrations.  Solution 
of  this  problem  is  accomplished  by  using  mathematical  .optimization 
techniques.  The  objective  function  is  related  to  meeting  downstream  target 
qualities  subject  to  various  hydraulic  constraints  on  the  individual  ports. 

Kaplan  (21]  solved  a  similar,  although  more  difficult,  problem  by 
including  in  the  constraint  set  upper  and  lower  bounds  on  the  relase 
concentration  of  each  water  quality  constituent.  Kaplan  also  considered  as 
part  of  his  objective  function  the  reservoir  water  quality  that  resulted  from 
any  particular  operation  strategy.  A  penalty  function  approach  was  used  to 
incorporate  the  many  constraints  into  the  objective  function  which  cuuid  then 
be  3olved  as  an  unconstrained  nonlinear  problem.  For  the  problem  ot  interest 
with  respect  to  HEC-5Q,  with  appropriate  transformations  it  is  possible  to 
formulate  a  quadratic  objective  function  with  linear  constraints. 
Mathematical  optimization  techniques  are  available  to  exploit  the  special 
structure  of  this  problem  and  to  solve  it  efficiently. 

The  hydraulic  structure  under  consideration  is  composed  of  two  wet 
wells,  containing  up  to  eight  ports  each,  and  a  flood  control  outlet.  It  is 
assumed  that  releases  through  any  of  these  ports  (including  the  flood  control 
outlet)  leave  the  reservoir  through  a  common  pipe.  At  any  given  time,  only 
one  port  in  either  wet  well  and  the  flood  control  outlet  may  be  operated. 
Hence,  the  algorithm  provides  flows  through  three  ports  at  most. 

The  HEC-5Q  model  also  provides  for  releases  through  an  uncontrolled 
spillway.  These  releases  are  not  a  part  of  the  gate  selection  algorithm,  but 
the  water  quality  of  the  spillway  releases  are  considered  by  the  gate 
selection  algorithm. 
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The  algorithm  proceeds  by  considering  a  sequence  of  problems,  each 
representing  a  different  combination  of  open  ports.  For  each  combination  the 


optimal  allocation  of  total  flow  to 
water  quality  index  is  determined  for 
combination  of  open  ports  with  the 
associated  allocation  of  flows,  define 
time  period  under  consideration. 


ports  is  first  determined  and  then  a 
the  optimal  allocation  of  flows.  The 
highest  water  quality  index  and  its 
the  optimal  operation  strategy  for  the 


There  are  four  different  types  of  combinations  of  open  ports.  For 
one-port  problems,  all  of  the  flow  is  taken  from  a  single  port  and  the  water 
quality  index  is  computed.  For  two-port  problems,  combinations  of  one  port 
in  each  wet  well  and  combinations  of  each  port  with  the  flood  gate  are 
considered.  For  three-port  problems,  combinations  of  one  port  in  eacn  wet 
well  and  the  floodgate  are  considered.  The  total  flow  to  be  released 
downstream  is  specified  external  to  the  port  selection  module,  but  if  the 
flow  alteration  option  is  selected,  then  the  flow  can  be  treated  as  an 
additional  uecision  variable  and  the  fl  w  for  which  the  water  quality  index 
is  maximized  is  also  determineu. 


For  each  combination  of  open  ports,  a  sequence  of  flow  allocation 
strategies  is  generated  using  a  gradient  method,  a  gradient  projection 
method,  or  a  Newton  projection  method  as  appropriate.  The  value  of  any  flow 
allocation  strategy  is  determined  by  evaluation  of  a  water  quality  index 
subject  to  the  hydraulic  constraints  of  the  system.  The  sequence  converges 
to  the  optimal  allocation  strategy  for  the  particular  combination  of  open 
ports. 


To  evaluate  the  water  quality  index  for  a 
strategy,  first  the  release  concentration  for  every 
is  computed. 


feasible  flow 
water  quality 


allocation 
const i tuent 


p 

p  =  1 

■i  0  1 

L  °!)  pi 

M 

‘  P 

i;  q 

P  -  1  11 

where : 

c  =  index  for  constituents 
Rc  =  release  concentration  for  constituent  c 
p  =  index  for  open  ports 
Np  ■=  number  of  open  ports 
***  . .  5  concent  ra  t  ion  of  constituent  c  at  port  p 
Qp  '  flow  through  port  p 

Nc  -  number  ot  constituents  under  consideration 
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The  deviation  of  release  qualities  from  downstream  target  qualities  can 
be  computed. 


D 

c 


R 

c 


;  c 


!.  Nc 


(221 


where : 

Dc  =  deviation  of  constituent c 

tc  =  downstream  target  quality  for  constituent c 

The  subindex  Sc  for  each  constituent  can  be  determined  bv: 

Sc  =  f  ( Dc )  ;  c  =  1,  Nc  (23) 


Where  the  function  f(Dc)  takes  the  quadratic  form: 

2 

f(Dc)  =  a  +  cDc  (24) 

Suggested  coefficients  for  some  water  quality  constituents  ar-'  shown  in 
Table  1. 

Table  1 

Coefficients  in  Constituent  Suboptimization 
of  Gate  Selection  Procedure 


Quality 

a 

c 

Temperature 

100 

-4.0 

DO 

100 

-4.0 

TDS 

100 

-0.00062b 

BOD 

100 

-0.444 

E.  Coli 

100 

-0.025 

NH3 

100 

-0.16 

N03 

100 

-0.16 

Finally,  the  scalar  water  quality  index  can  be  determined  by: 


N 

c 

z  =  >:  w  s 

-  c  c 

C=1 


( 2  S ) 
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wnere : 


Z  -  water  quality  index 

Wc  =  weigh; ing  factor  tor  constituent  e;  Lhe  sum  ul  the  weighting 
factors  for  all  constituents  must  equal  one 
Sc  =  ,ubindex  for  constituent  r 


In  summary,  the  problem  of  determining  the  optimal  allocation  ol  flows  to 
ports  tor  a  particular  combination  of  open  ports  and  for  a  specified  total 
flow  rate  Q,  can  he  expressed  as  follows: 


MAX 


o 

1’ 


c  =  l 


(26) 


Subject  to: 


"P 


.1  =  0 

V 


(27) 


m  i  n ,  p  =  Q 

P 


v 

'  max ,  p 


US) 


Wliere  i,'rnin  -stid  Fmax  ire  tne  minimum  and  maximum  acceptable  flow  rates 
through  a  port , 

When  an  acceptable  flow  range  Qj  ,wer  to  Qupper  ls  specified,  then  tne 
problem  is  written  as: 


MAX 

n 

'P 


W  s 


(29) 
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Subject  to: 


N 

0,  1  7.v  n  <  o 

lower  ,  p  upper 

i)  =  l 


(  30/ 


'  ,  '  n  i  7  ;  l'  = 

mln.p  |i  ni;i 


(31) 


These  problems  are  solved  very  efficiently  by  using  mathematical 
optimization  techniques  that  cake  advantage  of  tne  problem  structure,  name  1 v 
a  quadratic  objective  function  with  linear  constraints. 

Flow  Alterations 


Tne  flow  alteration  routine  is  designed  to  change  the  reservoir 
releases,  computed  by  the  flow  simulation  module,  to  better  satisfy  tne 
stream  control  point  water  quality  objectives.  The  routine  is  a<  signed  about 
a  mass  balance  for  all  reservoir  releases  and  all  control  points  afiected  Dy 
those  releases.  Tributary  inflows  and  other  flow  changes  ire  included. 
Second  order  effects,  such  as  reaeration  and  external  heating  due  to 
increased  or  decreased  stream  surface  area  are  not  included. 

The  procedure  is  at,  follows: 

1.  The  relative  mass  that  needs  to  be  added  in  the  flow  at  the  control 
point  (for  those  constituents  below  the  target)  or  reduced  in  the 
flow  at  the  control  point  (for  those  constituents  above  tne  target) 
is  computed  using: 


AM  =  Q  (C  -  C  )  (  32  ) 

cp  o  cp 


where : 

QCp  =  flow  at  the  controi  point  as  determined  by  the 
flow  simulation  module 

CQ  =  target  constituent  concentration  at  the  control 
point 

CCp  -  computed  constituent  concentration  at  the  control 
point 

2.  The  average  reservoir  release  concentration  is  computed  for  all 
reservoirs  for  which  the  constituent  concentration  in  the  releases 
is  greater  than  the  target  concentration  at  tne  control  point  ol 
interest  (for  those  constituents  below  the  target)  or  for  which  the 
constituent  concentration  in  the  releases  is  less  than  the  target  at 
the  control  point  of  interest  (for  those  constituents  above  the 
target ) .  Th  us : 


n,  .c.)/;  o 

.  .  1‘  i  i  .  .  -Ri 

i-l  1=1 


where : 

<Jg  -  average  constituent  concentrat ion  in  reservoir  releases 
for  only  those  reservoirs  releasing  flow  with 
constituent  concentrations  adequate  to  dilute  the 
control  point  concentration  and  bring  it  to  the  target 
QKl  =  flow  release  from  reservoir  i 

=  constituent  concentration  in  release  from  reservoir  i 
n  =  number  of  reservoirs 

and  the  sums  are  taken  only  over  those  reservoirs  i  that  are  capable 
ol  diluting  the  controL  point  constituent  concentration  that  is 
worse  than  the  target. 

The  total  dilution  flow  requirement  is  then  computed  by  the 
following  quotient: 


where  is  the  total  flow  release  needed  to  bring  the  ccnsti  cent 

concentration  at  the  control  point  of  interest  to  the  target. 

4.  The  flow  is  then  apportioned  to  the  reservoirs  capable  of 

ormging  the  control  point  constituent  concentration  to  the  target 
irs  proportion  to  the  flows  originally  computed  for  those  reservoirs 
by  the  flow  simulation  module. 

Tnus  the  flow  augmentation  requirement  can  be  computed  for  each  control 
point  and  for  each  constituent.  The  various  computed  flow  rates  are  then 
combined  by  using  the  coefficients  of  the  linear  programming  objective 
function  and  the  deviation  of  the  respective  constituent  concentrations  from 
the  target  concentrations  at  each  respective  control  point  as  follows: 


rp  re 


1-1  j  1 


N  N 
cp  cc 

i-l  .i-i 


P  (C  -C  I 
1  :  !  i  i" 


o,  1 

'k  i]  I)  10 


where : 


'Jk 

N'cp 


i  I 


»  J 


flow  release  from  reservoir  k 

number  ot  control  points  aftecteij  by  :>.»th  res  -r vo  1  r 
number  ut  constituents 

linear  programming  objective  function  toel  In  n  i:  f  •• 
constituent  j  at  control  point  i 
computed  concentration  of  constituent  i  at  coin  ■! 
point  i 

Lciigtf  L  COilCcIlll  a  l  luu  ui  C  Oiib  L  1  l  ticill  1 


Once  the  is  determined,  using  equation  (ib),  tiie  '  !  simulation 

module  is  recalled  and  the  daily  computations  lor  flow  and  wat- :  u  Mlitv  >r- 
solved  again  for  tlie  final  results. 


SUMMARY 

HEC-5Q  is  capable  of  simulating  the  effects  of  the  operation  f  as  nun/ 

as  ten  reservoirs  and  the  streams  of  the  basin.  The  reservoirs  may  :i. 

operated  to  satisfy  a  number  of  objectives,  including  flood  control,  low  flow 
maintenance,  hydropower  production,  water  conservation  and  water  quality 
control.  The  water  quality  portion  of  tile  model  will  simulate  temperature 
and  seven  water  quality  constituents  including  an  option  for  dissolved 
oxygen.  The  model  will  internally  determine  the  water  quality  needed  from 
all  reservoir  releases  to  meet  specified  downstream  water  quality  objectives 

and  will  determine  the  gate  openings  in  each  reservoir  that  will  yield  the 

appropriate  reservoir  release  water  quality.  Should  it  be  necessary,  flows 
will  be  altered  to  ensure  that  downstream  water  quality  objectives  are  met. 
As  currently  formulated,  the  model  does  not  use  foresight  in  an  attempt  to 
ensure  that  a  "global"  optimum  solution  is  found  that  meets  water  quality 
objectives.  The  model  selects  the  "best"  solution  for  system-wide  reservoir 
operation  on  a  daily  basis. 
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Appendix  -  Notation 


lhe  following  symbols  are  used  in  this  paper: 

A  -  surface  area  of  fluid  element 

,A|  -  matrix  of  coefficients  describing  water  elements  characteristics 
a  =  empirical  elements 

b  -  vector  of  constants  from  initial  conditions 

C  =  simulated  water  quality  concentration 

C0  =  target  water  quality  concentration 

c  =  specified  heat  of  water 

c  =  vector  water  quality  concentration 

D  =  effective  diffusion  coefficient 

Dc  -  deviation  of  actual  water  quality  concentration  from  target 
value 

Um  =  mol -Cellar  diffusion  coefficient 
d  •-  flow  depth 
K  -  acceptable  flow  rate 
H  =  external  heat  source  or  sink 


'*br 

Hc 

He 

^n 
II  . 


II 


s  r 


long  wave  atmospheric  r filiation 

iefl-ctea  iong  wave  radiation 

water  surface  radiation 

heat  of  conduction 

heat  of  eviporation 

net  iieat  transfer 

short  wave  solar  radiation 

reflected  short  wave  radiation 


Appendix  -  Nocation  (Continued) 


I  =  light  energy  at  specified  depth 
1D  =  light  energy  at  water  surface 
K  =  heat  exchange  coefficient 
=  decay  rate 
K.2  =  reaeration  rate 

k  =  light  extinction  coefficient 
m  =  number  of  rows  in  matrix 

Ncc  =  number  of  water  quality  constituents 

NCp  =  number  of  control  points  effected  by  nore  than  one  reservoir 

Np  =  number  of  open  ports  on  discharge  structure 

n  =  number  of  columns  in  matrix 

0  =  dissolved  oxygen  concentration 

0*  =  saturation  dissolved  oxygen  concentration 

p  =  index  for  open  ports  on  discharge  structure 
P^j  =  linear  programming  objective  functions 
Q  =  interelement  flow 
Qd  =  downward  advective  flow 
Qi  =  lateral  inflow 

=  reservoir  release 
Q0  *  lateral  outflow 
Qu  =  upward  advective  flow 
Qw  =  flow  withdrawal 
Qx  =  flow  into  element 
Rc  =  release  concentration 


Appendix  -  Nota' ion  (Continued) 

T  =  water  temperature 

rc  -  water  temperature  correction  factor 
Te  =  equilibrium  temperature 
Ts  =  water  surface  temperature 
Tx  =  inflow  water  temperature 
t  -  t ime 
II  =  flow  velocity 
V  =  volume  of  fluid  element 
Sc  =  water  quality  subindex 
W  =  wind  speed 

Wc  =  water  quality  index  weighting  factor 
Z  =  water  quality  index 

z  =  space  coordinate  (vertical  for  reservoir/hori zonta 1  for  stream) 
&z  =  element  thickness 
p  =  water  density 
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TP-56  Experiences  of  the  Hydrologic  Engineering 
Center  in  Maintaining  Widely  Used 
Hydrologic  and  Water  Resource 
Computer  Models,  B.  S.  Ei chert, 

Nov  78,  16  pp. 

TP-57  Flood  Damage  Assessments  Using  Spatial 

Data  Management  Techniques,  D.  W.  Davis 
and  R.  P.  Webb,  May  78,  27  pp. 


ADA-106  255 


ADA-106  256 


TP-58  A  Model  for  Evaluating  Runoff-Quality  in  ADA-  106  257 

Metropolitan  Master  Planning, 

L.  A.  Roesner,  H.  M.  Nichandros, 

R.  P.  Shubinski,  A.  D.  Feldman, 

J.  W.  Abbott,  and  A.  0.  Friedland, 

Apr  72,  81  pp. 


HEC 

NUMBER 

TP-59 

TP-60 

TP-61 

TP- 62 

TP-63 

TP-64 

TP- 6  5 

TP-66 

TP-67 

TP-68 


TECHNICAL  PAPERS  (TP ) (Cont inued ) 


TITLE 


HEC  NTIS 

PRICE  NUMBER 


$2.00  Each 


Testing  of  Several  Runoff  Models  on  an 
Urban  Watershed,  J.  Abbott, 

Oct  78,  53  pp. 

Operational  Simulation  of  a  Reservoir 
System  with  Pumped  Storage, 

G.  F.  McMahon,  V.  R.  Bonner  and 
B.  S.  Eichert,  Feb  79,  32  pp. 

Technical  Factors  in  Small  Hydropower 
Planning,  D.  W.  Davis,  Feb  79, 

35  pp. 

Flood  Hydrograph  and  Peak.  Flow  Frequency 

Analysis,  A.  D.  Feldman,  Mar  79  21  pp. 

HEC  Contribution  to  Reservoir  System 
Operation,  B.  S.  Eichert  and 
V.  R.  Bonner,  Aug  79,  28  pp. 

Determining  Peak-Discharge  Frequencies  in 

an  Urbanizing  Watershed:  A  Case  Study, 

S.  F.  Daly  and  J.  C.  Peters,  Jul  79,  15  pp. 

Feasibility  Analysis  in  Small  Hydropower 

Planning,  D.  W.  Davis  and  B.  W.  Smith, 

Aug  79,  20  pp. 

Reservoir  Storage  Determination  by  Computer 
Simulation  of  Flood  Control  and 
Conservation  Systems,  B.  S.  Eichert, 

Oct  79,  10  pp. 

Hydrologic  Land  Use  Classification  Using 
LANDSAT,  R.  J.  Cermak,  A.  D.  Feldman 
and  R.  P.  Webb,  Oct  79,  26  pp. 

Interactive  Nonstructural  Flood-Control 
Plannng,  D.  T.  Ford,  Jun  80,  12  pp. 


ADA- 106  258 

ADA-106  259 

ADA- 109  757 

ADA-109  758 

ADA-109  759 

ADA-109  760 

ADA-109  761 

ADA-109  762 

ADA-109  763 

ADA-109  764 
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HEC 

NUMBER 

TP-69 

TP-70 

TP-71 

TP- 7  2 

TP- 7  3 

TP-74 

TP-75 

TP- 76 

TP- 7  7 

TP-78 


TECHNICAL  PAPERS  ( TP) ( Cont inued ) 


TITLE 


HEC  NTIS 

PRICE  NUMBER 

$2.00  Each 


Critical  Water  Surface  by  Minimum  Specific 
Energy  Using  the  Parabolic  Method, 

B.  S.  Eichert,  1969,  15  pp. 

Corps  of  Engineers  Experience  with 
Automatic  Calibration  of  a 
Precipitation-Runoff  Model,  D.  T.  Ford, 

E.  C.  Morris,  and  A.  D.  Feldman, 

May  80,  12  pp . 

Determination  of  Land  Use  from  Satellite 

Imagery  for  Input  to  Hydrologic  Models, 

P. .  P.  Webb,  R.  Ceraak,  and  A.  D.  Feldman, 
Apr  80,  18  pp. 

Application  of  the  Finite  Element  Method  to 
Vertically  Stratified  Hydrodynamic  Flow 
and  Water  Quality,  R.  C.  HacArthur  and 
W.  R.  Norton,  May  80,  1?  pp . 

Flood  Mitigation  Planning  Using  HEC-SAM, 

D.  W.  Davis,  Jun  80,  17  pp. 


ADA-951  599 


ADA-109  765 


ADA- 109  766 


ADA- 109  767 


ADA-109  756 


Hydrographs  by  Single  Linear  Reservoir  ADA-109  768 

Model,  J.  T.  Pederson,  J.  C.  Peters, 
and  0.  J.  Helweg,  May  80,  17  pp. 


HEC  Activities  in  Reservoir  Analysis,  ADA-109  769 

V.  R.  Bonner,  Jun  80,  10  pp. 

Institutional  Support  of  Water  Resource  ADA-109  770 

Models,  J.  C  Peters,  May  80,  23  pp . 

Investigation  of  Soil  Conservation  Service  ADA-109  771 

Urban  Hydrology  Techniques, 

D.  G.  Altman,  W.  H.  Espey,  Jr.  and 
A.  D.  Feldman,  May  80,  14  pp . 

Potential  for  Increasing  the  Output  of  ADA-109  772 

Existing  Hydroelectric  Plants, 

D.  W.  Davis  and  J.  J.  Buckley, 

Jun  81,  20  pp. 
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TECHNICAL  PAPERS  (TP) (Continued ) 


HEC 

NUMBER 


HEC 

PRICE 


NTIS 

NUMBER 


P2.00  Each 


TP-79 


TP-80 


TP- 81 


TP-82 


TP-83 


TP-84 


TP- 8  5 


TP-86 


TP-87 


TP-88 


TP-89 


Potential  Energy  and  Capacity  Gains  from 
Flood  Control  Storage  Reallocation 
at  Existing  U.  S.  Hydropower 
Reservoirs,  B.  S.  Eichert  and 

V.  R.  Bonner,  Jun  81,  18  pp . 

Use  of  Non-Sequential  Techniques  in  the 

Analysis  of  Power  Potential  at  Storage 
Projects,  G.  M.  Franc,  Jun  81,  18  pp. 

Data  Management  Systems  for  Water  Resources 
Planning,  D.  W.  Davis,  Aug  81,  12  pp. 

The  New  HEC-1  Flood  Hydrograph  Package,  A.  D. 
Feldman,  P.  B.  Ely  and  D.  M.  Goldman, 
May  81,  28  pp. 

River  and  Reservoir  Systems  Water  Quality 
Modeling  Capability,  R.  G.  Willey, 

Apr  82,  15  pp. 

Generalized  Real-Time  Flood  Control  System 
Model,  B.  S.  Eichert  and  A.  F.  Pabst, 

Apr  82,  18  pp. 

Operation  Policy  Analysis:  Sam  Rayburn 
Reservoir,  D.  T.  Ford,  R.  Garland 
and  C.  Sullivan,  Oct  81,  16  pp. 

Training  the  Practitioner:  The  Hydrologic 
Engineering  Center  Program, 

W.  K.  Johnson,  Oct  81,  20  pp. 

Documentation  Needs  for  Water  Resources 

Models,  W.  K.  Johnson,  Aug  82,  16  pp. 

Reservoir  System  Regulation  for  Water 
Quality  Control,  R.G.  Willey, 

Mar  83,  18  pp. 

A  Software  System  to  Aid  in  Making  Real-Time 
Water  Control  Decisions,  A.  F.  Pabst 
and  J.  C.  Peters,  Sep  83,  17  pp. 


ADA- 109  787 


ADA- 109  788 


ADA-114  650 


ADA-114  360 


ADA-114  192 


ADA-114  359 


ADA- 123  526 


ADA-123  568 


ADA-123  558 


ADA-130  829 


HEC 

NUMBER 

TP-90 

TP- 91 

TP-92 

TP- 9  3 

TP-94 

TP-95 

TP- 96 

TP-97 

TP-98 

TP-99 


TECHNICAL  PAPERS  ( TP ) < Cont inued ) 


TITLE 


HEC 

PRICE 


$2.00  Each 


Calibration,  Verification  and  Application 
of  a  Two-Dimensional  Flow  Model, 

D.  M.  Gee,  Sep  83,  6  pp. 

HEC  Software  Development  and  Support, 

B.  S.  Eichert,  Nov  83,  12  pp. 

Hydrologic  Engineering  Center 
Planning  Models 
D.  T.  Ford  and  D.  W.  Davis, 

Dec  83,  17  pp. 

Flood  Routing  Through  A  Flat,  Complex 

Floodplain  Using  A  One-Dimensional. 
Unsteady  Flow  Computer  Program, 

J.  C.  Peters,  Dec  83,  8  pp. 


Dredged-Material  Disposal  Management 
Model,  D.  T.  Ford,  Jan  84,  18  pp . 


Infiltration  and  Soil  Moisture  Redistribution 
in  HEC- 1 ,  A.  D.  Feldman,  Jan  84, 


The  Hydrologic  Engineering  Center  Experience 
in  Nonstructural  Planning,  W.  K.  Johnson 
and  D.  W.  Davis,  Feb  84,  7  pp. 


Prediction  of  the  Effects  of  a  Flood  Control 
Project  on  a  Meandering  Stream, 

D.  M.  Gee,  Mar  84,  12  pp . 


Evolution  in  Computer  Programs  Causes  Evolution 

in  Training  Needs:  The  Hydrologic  Engineering 
Center  Experience,  V.  R.  Bonner,  Jul  84,  20  pp . 


Reservoir  System  Analysis  for  Water  Quality, 

J.  H.  Duke,  D.  J.  Smith  and  R.  G.  Willey, 
Aug  84,  2.  pp. 


NTIS 

NUMBER 


5 


